Phase-space analysis of a class of k-essence cosmology 
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We perform a detailed phase-space analysis of a class of k-essence cosmology. We find the critical 
points can be divided into three classes: points unstable but the model stable, both points and the 
model stable, points stable but the model unstable. Like the case of points unstable but the model 
stable, the case of points stable but the model unstable is not relevant from a cosmological point of 
view, though they can be late-time attractors for the universe. So in order to study the possible final 
state of the universe, it is important to investigate not only the stability of the critical points but 
also the stability of the model. The case of both points and the model stable can further be divided 
into two classes: points only presenting decelerating phases and points at which all decelerating, 
constant-speed, and accelerating phases can appear; the final state of the universe dependents on 
the potential. 
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CN ; I. INTRODUCTION 

O ! In the last decade a convergence of independent cosmological observations suggested that the Universe is experi- 
encing accelerated expansion. An unknown energy component, dubbed as dark energy, is proposed to explain this 
$^ \ acceleration. Dark energy almost equally distributes in the Universe, and its pressure is negative. The simplest and 
most theoretically appealing candidate of dark energy is the vacuum energy (or the cosmological constant A) with a 
constant equation of state (EoS) parameter w — —1. This scenario is in general agreement with the current astronom- 
er^ , ical observations, but has difficulties to reconcile the small observational value of dark energy density with estimates 
J> ' from quantum field theories; this is the cosmological constant problem P, [l]- Recently it was shown that ACDM 
model may also suffer an age problem 0] . It is thus natural to pursue alternative possibilities to explain the mystery 
of dark energy. Over the past decade numerous dark energy models have been proposed, such as quintessence, phan- 
tom, k-essence, tachyon, (generalized) Chaplygin gas, DGP, quintom [4], etc. k-essence, a simple approach toward 
constructing a model for an accelerated expansion of the Universe, is to work with the idea that the unknown dark 
^^O ■ energy component is due exclusively to a minimally coupled scalar field 4> with non-canonical kinetic energy which 
. results in the negative pressure [^j . 

K-essence scenario has received much attention, it was originally proposed as a model for inflation [6j , and then as 
a model for dark energy Q. In several cases, k-essence cannot be observationally distinguished from quintessence 
A method to obtain a phantom version of FRW k-essence cosmologies was devised in Q . The stability of k-essence 
k>( \ was studied in ^ . Dynamics of k-essence were discussed in [l3| . Conditions for stable tracker solutions for k-essence 
- • in a general cosmological background were derived in ll|. Slow- roll conditions for thawing k-essence were obtained 



' in [12l |. A connection between the holographic dark energy density and the kinetic k-essence energy density was 
discussed in fT3l|. An holographic k-essence model of dark energy was proposed in [l3]- The geometrical diagnostic 
for purely kinetic k-essence dark energy was discussed in jlS] . The equivalence of a barotropic perfect fluid with a 
k-essence scalar field was considered in [l^, [13] . Dark matter to dark energy transition in k-essence cosmologies was 
examined in [Tsll . Models of dark energy with p urely kinetic multiple k-essence sources that allow for the crossing of 
the phantom divide line were investigated in |19l |20| . The thermodynamic properties of of k-essence was discussed in 
[2l| . Models of k-essence unified dark matter were discussed in [22H26j . Power- law expansion in k-essence cosmology 
was investigated in 1271. Theoretical and observational Constraints on k-essence dark energy models were discussed 
in psi - lsoj . In Ref. [3 if, a model independent method of reconstructing the Lagrangian for the k-essence field by 
using three parametrizations for the Hubble parameter H{z) was studied in detail. With assumptions on the EoS of 
k-essence as functions of the scale factor a, Ref. )32] discussed the forms of the Lagrangians. By assuming the EoS 
of k-essence as functions of the kinetic energy, Ref. [3j] considered restrictions on purely kinetic k-essence. Ermakov 
invariant of k-essence and its observational probes for the early universe are discussed in [3^ . 
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Since from it was proposed, k-essence was been studied intensively. It is still worth investigating in a systematic way 
the possible cosmological behavior of the k-essence cosmology. Thus, in the present work we perform a phase-space 
and stability analysis of a class of k-essence cosmology. We are interesting in investigating the possible late-time 
solutions. In these solutions we calculate various observable quantities, such are the density of the dark energy and 
the EoS parameters. As we see, indeed k-essence cosmology can be consistent with observations. 

This paper is organized as follows, in the following section, we review the model of k-essence. In Sec. Ill, we 
perform a phase-space and stability analysis of k-essence cosmology. Finally, we shall close with a few concluding 
remarks in Sec. IV. 



II. K-ESSENCE COSMOLOGY 



Scalar fields with non-canonical kinetic terms often appear in particle physics. In general the action for such theories 
can be expressed as UHl 



S — d'^X\/—g 



(1) 



where p{(t),X) is a function in terms of a scalar field and its kinetic energy X = —^d^4)d^(f). We assume a flat 
and homogeneous Friedmann- Robertson- Walker (FRW) space-time and work in units SttG ^ c ~ 1. The action 
ffl includes a wide variety of theories, such as low energy effective string theory 6], ghost condensate 'sg'I, tachyon 
[37„i38,J, and DBI theory [33]. Here we consider an action which was extensively studied in literatures ^, i23„40 43] 



= V{^) [-X + X^] = V{^)F{X). 

where is the k-essence potential. For a constant potential Vo, action ([2]) are discussed in [23 

density of k-essence take the forms 

= Vm2XFx~F]^V{?>X^ -X), 
here Fx = dF/dX. The corresponding EoS parameter and the effective sound speed are given by 



(2) 

El]. The energy 



(3) 



2XFx - F 
2 dp/dX _ 



X'-X 
3X^~X' 
Fx 



_2X -I 
dp/dX ^ Fx + 2XFxx ^ 6X - V 



(4) 
(5) 



with Fxx = d^F/dX^. The definition of the sound speed comes from the equation describing the evolution of 
linear adiabatic perturbations in a k-essence dominated Universe [35| (the non-adiabatic perturbation was discussed 
in O, [4^ , here we only consider the case of adiabatic perturbations) . Perturbations can become classic unstable if 
the sound speed is imaginary, Cg < 0, so we insist on Cg > 0. Another potentially interesting requirement to consider 
is Cg < 1, which says that the sound speed should not exceed the speed of light, which suggests violation of causality. 



(6) 



Note, however, this is still an open problem (see e. g. |46l - l5l| ). 

In a flat and homogeneous FRW space-time, the equation for the k-essence field takes the form 

{Fx + 2XFxx)4> + iHFx^ + {2XFx - i^)-^ = 0, 



where = dV/ d(j), H is the Hubble parameter related to the Friedmann equations 



H^^^{Pm+P^), (7) 

H = -^{pm + P4>+P4>)- (8) 

Although we could straightforwardly include baryonic matter and radiation in the model, we neglect them for sim- 
plicity. In next section, we will transform Eqs. ([7]) and (jS]) into an autonomous dynamical system to perform the 
phase-space and stability analysis. 
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III. PHASE-SPACE ANALYSIS 



In general, in order to perform the phase-space and stabihty analysis, we have to transform the cosmological 
equations into an autonomous dynamical system by introducing the auxiliary variable x and y 0, [53 - l62j . Using the 
auxiliary variables x and y, the cosmological equations of motion, give a self-autonomous system X' = f(X), where X 
is the column vector constituted by the auxiliary variables, f(X) the corresponding column vector of the autonomous 
equations, and the prime denotes a derivative with respect to the logarithm of the scale factor, iV = Ina. The critical 
points Xc are extracted satisfying X =0, and in order to determine the stability properties of these critical points 
we expand around Xc, setting X = Xc + U with U the perturbations of the variables considered as column vector. 
Thus, up to the first order we acquire U' = M • U, where the matrix M contains the coefficients of the perturbation 
equations. Thus, for each critical point, the eigenvalues of M determine its type and stability. The condition for the 
stability of the critical points is: Tr M < and det M > 0. 



A. The basic equations and the critical points 



In order to transform the cosmological equations into an autonomous dynamical system, we have to introduce the 
auxiliary variable: 



(9) 



together with N = In a. It is easy to have — H-^. Using these variables we can straightforwardly obtain the 
density parameters of dark energy as 



^a:V(3x^-2). 



Because < fi^ < 1, we get the constraint equation 

1 



< ^x^y^iix^ - 2) < 1. 



(10) 



(11) 



in which the auxiliary variable x and y are allowed, see Fig. ([T]). In general, when the auxiliary variable x and y can 
take infinite values, it is necessary to analyze the dynamics at infinity by using Poincare Projection method [63l[63|. 
However, in the case we consider here, it is obviously that when |a:| > 4, y rapidly runs close to the zero; when y > 4, 
X rapidly runs close to -^2/3; When x — y take arbitrary values. So it is not necessary to analyze the dynamics at 
infinity in the case we discussed. 

The EoS, the sound speed of k-essence, and the total EoS are, respectively 



X 



2x2 



3x4 _ ' 

x2-l 



3x2 - 1 ' 

P0 



1 



xY{^'-2), 



A set of Eqs. 



and 



P4, + Pm 4 

gives a self-autonomous system in terms of the auxiliary variable x and y: 
, _ 3x(x2 - 1) V3 Ayx2(3x2 - 2) 



3x2-1 4 3x2-1 

l + \x^y\x^-2) 



(12) 

(13) 
(14) 

(15) 



y = — Y 2^ 



where the prime denotes a derivative with respect to the logarithm of the scale factor, N = Ina, and A = ~V^/V^ . 
Here we only consider the case in which the A is a constant, that is to say V{4)) cx (/)~2. Eq. (15) form an autonomous 
dynamical system in this case, this self-autonomous system are valid in the whole phase-space, not only at the critical 
points. 



4 



6 



5 

4 ^ 





FIG. 1: The allowed regions of the auxiliary variable x and y. 



The critical points {xc, Uc) of the autonomous system (|15p are obtained by setting the left hand sides of the equations 
to zero, namely let x' ~ y' ~ 0. After some algebraic calculus, we obtain the critical points as following: 
Critical point P\ 

xi = 0, (16) 
yi = 0. 

We calculate the values of Cg , fJ^, w^, and i^t, as shown in tablelH See for example: fi^ — 0, the universe is completely 
dominated by dark matter. We also note that in this case, remains unspecified and the results hold independently 
of its value. 

Critical point P2 

X2 = 1, (17) 
2/2 = 0. 

The corresponding cosmological parameters are presented in table ID Though in this case the k-essence behaves like 
the cosmological constant {w^ = —1), the universe is completely dominated by dark matter (fi^ = 0). 
Critical point P3 

X3 = -1, (18) 
2/3 = 0. 

In this case, the behavior of the universe is like that at point P2. Points Pi, P2, and P3 can be obtained in other 
scalar dark energy models (see e. g. 0), but the corresponding cosmological parameters, such as Cg , fi^, w^, and Wt, 
may be different because they are model-dependent. 
Critical point P4 

X4 = -\/2, (19) 
\/6 1 

with Wcf, = 0, meaning k-essence behaves like dark matter; £7^ = 3/A^ (A < — \/3), a universe is partly dominated by 
k-essence. 

Critical point P5 

^5 = %/2, (20) 
76 1 

= —A- 
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FIG. 2: Phase-space for k-essence cosmology, with the choice A = —3 for critical point Pi and A = 3 for critical point P5. 
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FIG. 3: Phase-space for k-essence cosmology, with the choice A = —1 for critical point Pe and A = 1 for critical point P7. 



In this case, the behavior of the universe for A > a/3 is like that at point P4. We plot critical points P4 for A = —3 
and P5 for A = 3 in Fig. H 
Critical point Pg 



/ A2 V48A2 + 9A4 

= -r+Y+ — 24 — 



(21) 



2/6 



(-20\/3A2 - 3V^A4 + 12V16A2 + SA* + 3AV16A2 + 3A4) 



4A(6 + A2) 



with = 1, meaning a k-essence dominated universe; and Wt = — 1 — + V48A^+9Af ^ (— VS < A < 0), all 
decelerating, constant-speed, and accelerating phases can appear. 
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FIG. 4: Phase-space for k-essence cosmology, with the choice A = — 1 for critical point Ps and A = 1 for critical point Pg. 



Critical point 



/ A2 V48A2 + 9A4 
-7 = \/l + ^ + , (22) 

X7(-20V3A^ - 3V3A4 + 12V16A2 + SA^ + 3AV16A2 + 
~ 4A(6 + A2) 

In this case, the behavior of the universe for < A < \/3 is like that at point Pq. We plot critical points Pe for A = — 1 
and P^ for A = 1 in Fig. H 
Critical point Pg 



v/l44 + 18A2 - 3V32A2 + W 
x& ^ ^ . (23) 



(-20V3A2 - 3V3A4 - 12V16A2 + 3A4 - 3AV16A2 + 3A4) 



4A(6 + A2) 

with 17^ = 1, a k-essence dominated universe; Wt — — 1 — ^A^ — V^sa^+qa ^ ^ q-j^ accelerating phases can 
appear. 

Critical point Pg 



\/l44 -I- 18A2 - 3V32A2 + W , , 

a^g = , (24) 

xg{-2{)^/l\^ ~ 3V3A^ - 12V16A2 + SA^ - 3A2V16A2 + 3A4) 
~ 4A(6-I-A2) 

In this case, the behavior of the universe for A < is like that at point P%. We plot critical points Pg for A = 1 and 
Pg for A = -1 in Fig. [H 



B. Stability analysis and the attractors 



Now we analyze stability of the critical points obtained above. For hyperbolic critical points (all the eigenvalues 
have real parts different from zero) one can easily extract their type (source (unstable) for positive real parts, saddle 
for real parts of different sign and sink (stable) for negative real parts). However, if at least one eigenvalue has a 
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TABLE I: The cosmological parameters and the behaviors of the critical points. 



zero real part (non-hyperbolic critical point) one is not able to obtain conclusive information about the stability from 
linearization and needs to resort to other tools like Normal Forms calculations ,65> ,66j , or numerical experimentation. 
The 2x2 matrix M of the linearized perturbation equations writes: 



M 



o , /^\ ^ Xxi (3x^-2) 



According to the conditions for the stability of the critical points: Tr M < and det M > 0, we obtain the ranges 
of A to make the critical points obtained above stable, as shown in table H] in which we also present the necessary 
conditions for their existence, as well as the corresponding cosmological parameters, c^, fl^, ly^, and Wt, with these 
parameters we can determine the stability of the model and discuss whether there exists acceleration phase or not. 

It is noteworthy that the stability of points do not mean the model is stable. So we must also investigate the 
stability of the model- both classical and quantum. The classical fluctuations may be regarded to be stable when the 
speed of sound is positive, Cg > 0. Combining this condition and Eq. (jl3p . we obtain the range in which the model is 

classic stable: X>i(x>lora::< -1), ox X < \ (-^ < a: < ^). 

In order to discuss the quantum stability of the scalar field, we consider the small fluctuations (50 around a back- 
ground value 00 which is a solution of the equations of motions: = 0o + ^4>- Since we are most concerned with 
UV instabilities, it is not restrictive to consider a Minkowski background metric and to choose, at least locally, a 
time direction in such a way satisfy the decomposition. By expanding p at the second order in (5(/), we can find the 
Lagrangian and the Hamiltonian fluctuations |42| : 

dU^ [px +^Xpxx)^ \-px ^ , (25) 

where = cPp/dcf)^. The positivity of the first two terms in Eq. (125^ leads to the following stability conditions 

px + 2Xpxx > 0, px> 0. (26) 

Because V{4>) > 0, we obtain X > 1/2 (x > 1 or x < —1) from these quantum stability conditions. So we conclude 
only for X > 1/2 (a; > 1 or x < —1) the model is both classic and quantum stable. With these conditions we can 
determine whether the model is stable or not at critical points when variations x and y take the corresponding values 
of critical points. 

In sum, the stability of points is related to the perturbations 6x and Sy, and depends on the condition: Tr M < 
and detM > 0; the classic stability of the model is related to the perturbations Sp, and depends on the condition: 
Cg > 0; and the quantum stability of the model is related to the perturbations 6(j), and depends on the condition: 
Px + '^Xpxx > and px > 0. That is to say, the stability of the points is different from the classic/quantum stability 
of the model. So it is important to investigate not only the stability of points but also the stability of the model. 
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From Table m it is obvious that points Pi, P2, and P3, are unstable for all A; points P4, Pq are stable for a 
certain range of A. The model is classic stable at points Pi, P7 for all A, classic unstable at points Pg and Pg for 
A 7^ 0; qua ntum stable at points P2, P7 for all A, quantum unstable at point Pi for all A, at points Ps and Pg for 
< A < VWl5- 

C. Cosmological implications 

According the phase-space analysis of a class of k-essence cosmology, we can now discuss the corresponding cos- 
mological behavior. The critical points Pi, P2, and P3, corresponding to complete dark matter dominated universe, 
are not relevant from a cosmological point of view, since they are unstable, though the model is classic and quantum 
stable at points P2 and P3 (quantum unstable at Pi). 

Points P4 with A < — -s/S and P5 with A > -s/S are more interesting because not only the points are stable but also 
the model are both classic and quantum stable. Thus they can be the late-time attractor of the universe. In this case, 
k-essence behaves like dark matter; the final state seems like a dark matter dominated universe. 

The critical points Pg with — \/3 < A < and P7 with < A < y/3 consist complete k-essence dominated and stable 
late-time solutions, both presenting accelerating phases if —^/3/3 < A < for Pg and < A < -\/3/3 for P7. The 
model is also both classic and quantum stable. For Pg with — \/3 < A < — \/3/3 and P7 with \/3/3 < A < -\/3, these 
two points present decelerating phases. For Pg with A = — \/3/3 and P7 with A = \/3/3, these two points result to 
constant-speed expansion. 

The critical points Pg with A > and Pg with A < also consist complete k-essence dominated and stable late-time 
solutions, both presenting accelerating phases, at which k-essence behaves like phantom. But the model is classic 
unstable and quantum unstable (0 < A < -y/lB/lS) at these points. So they are also not physical interesting. 

IV. CONCLUSIONS AND DISCUSSIONS 

In this work we have made a comprehensive phase-space analysis of a class of k-essence cosmology. We examined 
if a universe governed by k-essence can have late-time solutions compatible with observations. 

We obtained the critical points and presented the conditions for their existence and stability. We also calculated the 
values of the corresponding cosmological parameters, Cg, ft^, w^, and Wt, which are important in k-essence cosmology. 
We investigated the classic stability of the model, as well as the quantum stability. We discussed the behaviors of 
the critical points and plotted some of them. The critical points can be divided into three classes: unstable points at 
which the model is stable (for example, points P2 and P3), stable points at which the model is stable (for example, 
points P4, P5, Pg, and P7), stable points at which the model is unstable (classic or quantum unstable, for example, 
points Pi, Pg, and Pg). The case of points unstable but model stable or the case of points stable but model unstable 
are not relevant from a cosmological point of view. Only stable points at which the model is also stable (classic and 
quantum) are physically interesting. So only points P4, P5, Pg, and P7 are cosmological relevant. For points P4 and 
P5, which can be the late-time attractor for the universe, at which k-essence behaves like dark matter, and the final 
state seems like a dark matter dominated universe. For points Pg and P7, which can also be the late-time attractor for 
the universe at which k-essence completely dominates, but all decelerating, constant-speed, and accelerating phases 
can appear, which dependents on the potential. In this case, we still can not determine what is the final state of 
the universe. Although points Pg and Pg can indeed be the late-time attractor for the universe at which k-essence 
completely dominates and only accelerating phases appear, they are not physical concerning, because at which the 
model is both classic and quantum unstable. 

The results we obtained have shown that the potential of k-essence is important for the behavior of the universe. In 
order to study the possible final state of the universe, it is important to investigate not only the stability of the critical 
points but also the (classic and quantum) stability of the model. To our knowledge, no previous work considers both 
of them at the same time. The case combining the stability of the points and non-adiabatic perturbation of the model 
is physical interesting and will be consider anywhere. 

The analysis we performed indicates that a class of k-essence cosmology discussed here can be compatible with 
observations. But it just faces the problem from the cosmological point of view, and thus its results can been taken 
into account only if k-essence passes successfully observational tests which is the subject of interest of other studies. 
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